function u = CrankNicholson_diffusionConvection(u,N,x_max,k,p,num)
for i = 1:num
    u = CrankNicholson_diffusionConvection_onestep(u,N,x_max,k,p);
end